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It was recently shown that a black hole (BH) is the only compact object that can acquire a scalar 
charge in scalar Gauss-Bonnet (sGB) theory under the small coupling approximation. This leads to 
the fact that scalar radiation is emitted from a binary containing at least one BH. In this letter, we 
find the constraints on this theory from BH low-mass X-ray binaries (BH-LMXBs) . The main result 
of this letter is that from the orbital decay rate of A0620-00, we obtained a conservative bound that 
is six orders of magnitude stronger than the solar system bound. In addition to this, we look at 
XTE J1118+480, whose orbital decay rate has been recently measured with an excess compared to 
the theoretical prediction in GR due to the radiation reaction. The cause of this excess is currently 
unknown. Although it is likely that the cause is of astrophysical origin, here we investigate the 
possibility of explaining this excess with the additional scalar radiation in sGB theory. We find 
that there still remains a parameter range where the excess can be explained while also satisfying 
the constraint obtained from A0620-00. The interesting point is that for most of other alternative 
theories of gravity, it seems difficult to explain this excess with the additional radiation. This is 
because it would be difficult to evade the constraints from binary pulsars or they have already been 
constrained rather strongly from other observations such as solar system experiments. We propose 
several ways to determine whether the excess is caused by the scalar radiation in sGB gravity 
including future gravitational wave observations with space-borne interferometers, which can give a 
constraint three orders of magnitude stronger than that from A0620-00. 



Introduction: Testing gravitational theories [l| is im- 
portant from both theoretical and phenomenological 
points of view. For the former case, if we assume that 
the classical gravitational theory appears at the low- 
energy limit of a more fundamental theory such as su- 
perstring theory [3, H|, the theory does not necessarily 
reduce to general relativity (GR). One possibility of an 
effective gravitational theory is Einstcin-Dilaton-Gauss- 
Bonnet (EDGB) theory 0], where at the level of the 
action, the dilaton is coupled to the Gauss-Bonnet in- 
variant with a coupling constant a. For the latter case, 
modification of gravity may naturally solve problems in 
GR such as dark energy, dark matter and inflation. (See 
e.g. Ref. |5(.) So far, GR has been tested mainly in 
the solar system and in binary systems, especially binary 
pulsars [l[ . We have recently shown Q that only a black 
hole (BH) can acquire a scalar charge in scalar Gauss- 
Bonnet (sGB) gravity Q under the small coupling ap- 
proximation. Hence there exists scalar dipole radiation 
in a compact binary system where at least one of the 
constituents is a BH. (Here, sGB gravity refers to a the- 
ory where, at the level of the action, the Gauss-Bonnet 
invariant is coupled to an arbitrary function of a scalar 
field. EDGB theory is one specific theory in sGB gravity.) 
This means that we would not be able to test this theory 
with a neutron-star (NS)/NS or a white-dwarf (WD)/NS 
binary. In this letter, we probe this theory using BH low- 
mass X-ray binaries (BH-LMXBs). 
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One of the current bounds on sGB theory has been 
obtained using the Saturn probe Cassini Q, which mea- 
sured the Shapiro time delay, giving the constraint 
•v/H < 8.9 x 10 n cm. Especially, for EDGB theory, 
a stronger constraint has been obtained from the exis- 
tence of a stellar-mass BH Q. For example, the LMXB 
GRO J0422+32 is likely to contain a primary BH with 
mass (4 ± 1)M© [13] . Since a BH can only exist if \/\a\ 
is below the upper bound that is proportional to the 
Schwarzschild radius of a BH [ll[, we get the 2-er con- 
straint y/\a\ < 3.1 x 10 5 cm. 

LMXBs have been used to place constraints fl2l ITp 
on several theories such as Brans-Dicke (BD) theory UM . 
Especially, BH-LMXBs have been recently exploited (la . 
HH to put bounds on the size of the extra dimension in 
the Randall-Sundrum II braneworld model [llj ■ We first 
use the orbital decay rate of A0620-00 to obtain a new 
constraint on sGB gravity and find that it is more than 
six orders of magnitude stronger than the current bound 
from solar system experiments [13]. Among BH-LMXBs, 
XTE Jl 118+480 is an extremely interesting source since, 
very recently, Gonzalez et al. [l9[ reported that there is 
an excess in its orbital decay rate compared to that pre- 
dicted by GR. According to their arguments, this excess 
cannot be explained by the relativistic periastron preces- 
sion, rclativistic jets and the presence of a circumbinary 
disk. In order to interpret the cause as the magnetic 
braking, it is required that either the mass accretion 
onto the BH is almost zero (which is rather implausi- 
ble) or that the magnetic field of the secondary be more 
than 1-2 orders of magnitude greater than that typical in 
highly- rotating low-mass stars [13]. Although it is most 
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likely that this excess comes from an astrophysical origin, 
here, we investigate whether it can also be explained by 
the modification of gravity. It would be difficult for most 
alternative theories of gravity since they are already con- 
strained rather strongly, especially from binary pulsar ob- 
servations. However, sGB gravity cannot be constrained 
from these observations. Indeed, we find that the pos- 
sibility that the excess may be caused by the additional 
scalar radiation in sGB gravity has not been ruled out 
yet. We propose several ways to distinguish whether the 
excess is caused by this additional radiation. We use the 
geometrical unit c = G = 1 throughout the letter. 

sGB gravity: Let us consider the following action for 
sGB gravity 0: 



(1) 



Here, k = (167r) _1 , g stands for the determinant of the 
metric g^ v and R is the Ricci scalar. Rq B is the Gauss- 
Bonnet invariant defined as Rq B = R 2 — AR^R^" + 
R^upaR^ 1 " 7 with R pu and R^ vpa representing the Ricci 
and Riemann tensors, a is the coupling constant of the 
theory while £ ma t corresponds to the matter Lagrangian 
density. /(</>) is an arbitrary function of the scalar field 
<j> and f(<j>) = represents EDGB theory [S |20|. In 
c = G = 1 units, the scalar field <f> is dimcnsionlcss while 
the coupling constant a has a unit of (length) 2 . Following 
Ref. @ , we expand /(</>) around the asymptotic value <j> = 
at spatial infinity as f{<j>) = /(0) + f'(0)(j) + O((j) 2 ). The 
first term does not modify GR while /'(0) in the second 
term can be absorbed to a, hence we assume f(<j>) = <p. 

The action shown in Eq. (JTJ) should be treated as an 
effective theory: a truncation of some more fundamental 
or complete theory (such as supcrstring theory) to second 
order in a curvature expansion. It is only valid when 
the quadratic curvature term is smaller than the linear 
curvature term (i.e. a(/)i? 2 j B /(«J?) < 1). In order to 
obtain the small coupling condition that can be applied 
to a vacuum spacetimc as well, it would be convenient 
to extend this to acf>K 2 /(KK) = atfiK/n < 1, where K 2 
is defined as K 2 = G iivpa C iLVpa with C^ vpa representing 
the Weyl tensor. The correction relative to GR becomes 
the largest at the smallest length scale of the system. For 
a LMXB, this corresponds to the horizon radius of a BH, 



where (p ~ aj m 2 [2l| and K 



with mi representing 



the mass of the BH. Therefore, the condition that the 
effective theory remains to be valid is given by a<f>K j k ~ 
a 2 /(Kmf) < 1. Hence, throughout this letter, we apply 
the small coupling approximation where £ = a 2 / (nrnf) < 
1 . This is a reasonable approximation given that GR has 
passed many tests in the weak field regime. We only 
keep to linear order in £, which means that we only take 
4> up to ©(C 1 / 2 ). We can safely neglect the higher-order 
scalar field kinetic terms that appear in e.g. Ref. fl8| . 
The small-coupling approximation guarantees that the 



field equations to remain quadratic order in derivatives, 
which leads to the system being stable. 

In this theory, an object acquires a scalar charge pro- 
portional to the volume integral of the Gauss-Bonnet in- 
variant -Rq B @. Since Rqq is a topological invariant, 
the volume integral vanishes for any simply connected, 
asymptotically flat spacetimc. This means that the scalar 
charge vanishes for a star, including a NS, and only a BH 
can acquire a no n- vanishing scalar charge. To be precise, 
the scalar charge of a star does not exactly vanish since 
the universe is not simply connected due to other BHs. 
However, its effect can be safely neglected as long as its 
mass length scale is much smaller than the distance to 
its nearest BH. 

Orbital decay rate of a LMXB: Let us consider a cir- 
cular LMXB consisting of a BH with mass mi, a main 
sequence star with mass 777,2 and radius R2 and having a 
binary separation a. In Ref. |(|, we derived the energy 
flux E due to GW and scalar radiation, where the ori- 
gin of the leading correction relative to GR is the scalar 
dipole radiation and the order is of "-1 post-Newtonian 
(PN)". By using the relation E = VlL that holds under 
a circular-orbit binary, with Q = \Jmja? denoting the 
orbital angular velocity, we derive the decay rate of the 
orbital angular momentum L = jjL^Jma as [44| 



L = Lqjx (I + Av a ) . 



(2) 



Here, /1 is the reduced mass, m denotes the total mass, 
v = \J m/a is the typical velocity of the binary con- 
stituents, the parameters A and a are given by [45| 



A= — C. 
96^ 



a = -2. 



(3) 



and Lgr = —(32/5)ri 2 mv 7 with 7/ = niim^/m 2 repre- 
senting the symmetric mass ratio. 

Eqs. and ^ can be safely used for LMXBs when- 
ever the small coupling approximation is valid. 

The orbital period P = 2ir^/ 'a 3 Jm changes with time 
due to GW radiation and mass loss from the system by 
the stellar wind [46[ . The decay rate of the orbital period 
P for a LMXB can be estimated as QjL LUl 



P \D) L 



(4) 



where 77 and D are given in Ref. [13| in terms of the fol- 
lowing parameters: the mass ratio q = mi/m^, the mass 
transfer efficiency fi. t - n = —7771/777,2, the specific angular 
momentum j w carried away by the stellar wind in units 
of 2ira 2 /P and the adiabatic index £ a d = d In i?2/rfm 777,2. 
We assume that ftm and j w both take values from to 1 , 
while we fix £ a d as Cad =0.8 0,0. 

A new constraint on sGB gravity from A0620-00: Now, 
we derive the upper bound on £ from the system A0620- 
00 whose parameters are summarized in the first row of 
Table HI By using the observed bound on P and Eqs. @ 
and (13}, we obtain the 2-er upper bound on ( as ( < 
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TABLE I: The observed parameters (the primary mass mi, the mass ratio q, the orbital period P, the orbital decay rate P, 
the eccentricity e and the inclination i) of LMXBs A0620-00 and XTE J1118+480. 



System 


mi (M Q ) 


q 


P (hr) 


P (s/s) 


e 




A0620-00 
XTE J1118+480 


6.6 ±0.25 [22] 
8.30+°f 4 [19J 


17 ± 1 [15] 
37.0 ±12.3 [19] 


7.8 [15] 
4.08 [19] 


(1.66 ±2.64) x 10" 11 [15J 
-(5.8 ±2.1) x 10~ n [19] 


< 0.0067 [19] 


51° ±0.9° [22] 
68° ± 2° [231 



7.3 x 1CP 2 . (Here, we maximized £ over (i) /?m and j w 
in the range to 1, and (ii) P within 2-er range.) This 
leads to the 2-er constraint on y/\a\ as 

y/\a\ < 1.9 x 10 5 cm. (5) 

This new constraint is more than six orders of magnitude 
stronger than the solar system bound [lH. Eq. §5$ is 
the main result of this letter. Even if we restrict our 
attention to EDGB theory, it is still slightly stronger than 
the bound from GRO J0422±32. 

One might think that the constraint C < 7.3 x 1CP 2 
validates our assumption of small coupling approxima- 
tion only marginally. However, we cannot completely 
rule out the possibility that the small coupling approxi- 
mation is violated because at this stage, it would be diffi- 
cult to estimate the amount of contribution coming from 
the higher curvature terms that we have truncated. If 
the small coupling approximation is violated, it is cither 
the case where (i) higher curvature contribution is larger 
than 0(C) one, or that (ii) these two give comparable con- 
tributions. For the former, there would be much larger 
corrections relative to GR, and hence Eq. (|5|) still holds 
as a conservative bound. For the latter, if the signs of the 
two contributions are opposite, 0(C) correction would be 
reduced by the higher curvature one and the bound be- 
comes weaker than Eq. (|5|). However, notice that the 
bound on \/\a\ scales with A 1 / 4 (see Eqs. © and (J3|)). 
This means that even if the correction term in Eq. ([2]) 
is reduced by a factor of 10, the bound on i/|a[ is only 
weakened by a factor of 2. Moreover, this is a rather fine- 
tuned case compared to the case (i) or the one where the 
small coupling approximation is valid. For these reasons, 
we expect that the order of magnitude of the conserva- 
tive bound < 0(lO 5 )cm should still hold even if the 
small coupling approximation is not valid. 

Explaining the excess in the orbital decay rate of 
XTE J1118+480 in sGB gravity: For a LMXB XTE 
J1118±480 (whose observed parameters are summarized 
in the second row of Table H]) , the observed orbital de- 
cay rate is more than 10 times larger than the upper 
bound predicted in GR as |P| < 2.0 x 10~ 12 s/s. The 
cause of this excess is most likely to be of astrophysi- 
cal origin. However, we want to point out that this ex- 
cess can also be explained by the additional scalar radi- 
ation in sGB gravity. In order to explain this observed 
value, we found that -\/ja[ is required to be in the range 
1.1 x 10 5 cm < \/\q\ < 3.4 x 10 5 cm. (Here, we have both 



minimized and maximized \f\o\ over j w , P and q, 
where the range of the allowed value for q is shown in 
Table HI We have included q since its determination ac- 
curacy is not so high compared to other parameters.) 
Combining this with Eq. ([5]), we conclude that if \J\ot\ 
lies in the range 

1.1 x 10 5 cm < y/\c7\ < 1.9 x 10 5 cm, (6) 

the excess can be explained by the scalar radiation in 
sGB gravity. If we fix q = 37, the only difference in 
Eq. ([5]) is that the lower bound changes to 1.2 x 10 5 cm. 

The allowed range in Eq. ([6]) is rather small, but the 
point here is that at this stage, we cannot completely 
rule out the possibility of the excess due to the addi- 
tional scalar radiation. Below, we propose several ways 
to further pin down the value of a or that can completely 
rule out this possibility. 

Ways to test sGB gravity as an explanation for the ex- 
cess: It would be interesting to test Eq. ^ with future 
observations. One possibility is to look at the orbital de- 
cay rates of other BH-LMXBs since if the excess is caused 
by the scalar radiation, this effect should be universal to 
all similar systems. Alternatively, if the measurement ac- 
curacy of P in A0620-00 could be improved, it might be 
possible to probe the range found in this letter. How- 
ever, in order to realize this goal, the observation accu- 
racy must be improved by about an order of magnitude, 
which is quite challenging. In principle, if P can be mea- 
sured in the future, it may be possible to determine the 
cause of this excess (whether it is caused by the scalar 
radiation in sGB gravity, by the magnetic braking or by 
some other effect), since the correction in the energy flux 
should appear at different PN orders. However, this can 
only be done when the measurement accuracy of P im- 
proves considerably and this approach is much harder 
than the one previously mentioned. 

Another possibility is that future GW interferome- 
ters may be able to distinguish it. In Ref. Q, we 
found that the second-generation ground-based detectors 
such as adv. LIGO [24j are likely to place a constraint 
y/\c7\ < 4 x 10 5 cm with a BH binary of (6 ± 12)M Q and 
a signal-to-noise ratio (SNR) of 20. Unfortunately, this 
is not sufficient to probe the values shown in Eq. 
(and also the constraint is beyond the validity of the 
small coupling approximation) . For the Einstein Tele- 
scope (ET) [23, which is roughly 10 times more sen- 
sitive than adv. LIGO, the constraint is still slightly 
larger than the upper bound in Eq. ([5]) provided the 
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upper bound on y[a[ scales as (SNR) 1 / 4 . Next, let 
us consider space-borne GW interferometers. By per- 
formi ng ( sky-averaged) Fisher analyses |26l explained in 



Refs. 27, \l 



we found that ELISA [2J 
~ < 1.0 x 10 5 cm for a (10 + 



J can give the 
10 5 )M G 



constraint y/\a\ < l.U x lircm tor a (1U + LU")Mq circu- 
lar, spin-aligned BH binary at lGpc for lyr observation. 
Furthermore, DECIGO/BBO [HtH will be able to con- 
strain < 1.4 x 10 3 cm for a (1.4 + 10)M Q NS/BH 
binary with SNR = 10 for lyr observation. As discussed 
in Ref. [Hj|, DECIGO/BBO is expected to detect about 
N ~ 10 5 NS/BH binaries. Following Ref. H3], we expect 
that these multiple-source detection would make the con- 
straint of £ stronger by roughly y/~N. Since v/Jaf oc C 1 ^ 4 , 
we expect that the constraint on a can be improved to 
yja[ < 10 2 cm. This is indeed three orders of magni- 
tude stronger than the new constraint from A0620-00 
(Eq. ||5j|). We also expect that the results are almost the 
same for precessing binaries [13, [32| . If we include the 
eccentricity e, it would correlate strongly with £ since 
the correction terms from sGB and eccentricity both be- 
come greater for larger separation. However, for the same 
reason discussed in Ref. [23], we have checked that this 
degeneracy is expected to be solved by adding prior in- 
formation (e 2 > 0) to the eccentricity. 

Discussions: It seems difficult to explain the excess 
in the orbital decay rate of XTE J1118+480 in the 
context of most alternative theories of gravity, includ- 
ing massless and massive Brans-Dicke [3, [H| , dynami- 
cal Chern-Simons (CS) [34j . (non-linear) massive gravity 
theories l35l . |3(I , Einstein- Aether [53, HH and Hofava- 
Lifshitz [39|, |40| gravities, since most of these theories 
have already been constrained rather strongly from bi- 
nary pulsars and other tests such as solar system exper- 
iments [|| and galactic observations [4l|. In particular, 
even though it has not been constrained strongly yet, it 
would be difficult for the theories mentioned above to 
satisfy the requirements from the binary pulsar tests and 
also explain the excess in the LMXB. This is because the 
former observations are more precise than the latter and 
are known to agree with GR up to 0(1O -2 ). The interest- 
ing point about sGB gravity is that the scalar radiation 
appears only from binaries that contain BHs. This means 
that sGB gravity cannot be constrained from binary pul- 
sar observations. Hence, it can explain the excess within 



the current bound. 

In this letter, we have used the values of binary param- 
eters that have been determined in GR. This can be justi- 
fied as follows: For sGB gravity, in Ref. 0, we discussed 
that there should be two effects that give non-dissipative 
corrections to GR, (i) the scalar forces and (ii) the defor- 
mation in the metric. Since the secondary stars do not 
have scalar charges, we can forget about the former ef- 
fect. For the latter, there would be a 2PN correction, as 
found by Yunes and Stein [U, of O ((v 4 ) < 10" 11 , which 
is way below the measurement error in mass. Therefore, 
we can safely neglect the effect of modifications to the 
gravitational theory on the measured values of binary 
parameters. 

In Ref. @ , we also found the dissipative correction in a 
compact binary system in CS gravity. Therefore, it would 
be interesting to obtain new constraints on this theory 
from current binary observations. We expect that the 
double-pulsar binary [42j imposes a stronger constraint 
than LMXBs since the observation accuracy is higher for 
the former case. In order to investigate this goal, we first 
need to derive the conservative correction in this theory, 
which comes from the deformation in the metric at the 
quadratic order in spins and a dipole-dipole force that 
acts on each of the binary component. We expect that 
the double-pulsar binary constraint should be stronger 
than the current solar system bound [43j . This work is 
currently in progress. 

It would be interesting to extend our formalism de- 
veloped in Ref. @ for eccentric binaries and for non- 
vanishing scalar potential, and apply them to various bi- 
nary systems. These issues are left for future work. 
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